Naval  Research  Laboratory 

Washington.  DC  20375-5000 


NRL/MR/4404-92-6961 


Computations  of  Detonation  Structure: 
The  Influence  of  Model  Input  Parameters 


Michel  H.  Lefebvre 

Department  of  Chemistry 
Polytechnic  Institute 
Royal  Military  Academy 
Brussels,  Belgium 

Elaine  S.  Oran  and  K.  Kailasanath 
Laboratory  for  Computational  Physics  and  Fluid  Dynamics 


April  10,  1992 


92-09908 

’aiiHi] 


Approved  for  public  release;  distribution  unlimited 

4  ail 


REPORT  DOCUMENTATION  PAGE 


form  Approved 
0MB  No  OTC't-OIBB 


Public  reponir»9  burden  lor  thu  coHeoion  of  informjtion  estimated  to  average  t  hour  per  response,  including  the  time  lor  reviewing  imtructions.  searching  eaistmg  data  sources, 
gathering  and  maintaining  the  data  needed,  arsd  complettrtg  ar>d  reviewing  the  collection  ol  information  Send  comments  regarding  this  burden  estimate  or  any  other  asoect  of  this 
collection  of  information.  irKludmg  suggestions  for  reducirsg  this  burden,  to  Washington  Headquarters  Services.  Directorate  for  information  Operatiorn  and  Reports,  121S  Jefferson 
Davis  Highway.  Suite  1204.  Arlington.  VA  22202*4302.  and  to  the  Office  of  Management  arid  ftudget.  Paperwork  Reduction  Project  (0704*0188),  Washington.  DC  20S03 


1.  AGENC'^  USE  ONLY  (Leave  biank) 


4.  TITLE  AND  SUBTITLE 


REPORT  DATE 

April  10,  1992 


3.  REPORT  TVPE  ANO  OATES  COVERED 
Interim 


Computations  of  Detonation  Structure:  The  Influence  of  Model 
Input  Parameters 


6.  AUTHOR(S) 

Michel  H.  Lefebvre,  Elaine  S.  Oran,  and  K.  Kailasanath 


7.  PERFORMING  ORGANIZATION  NAME(S)  ANO  AOORESS(ES) 


Naval  Research  Laboratory 
Washington,  DC  20375-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/MR/4404-92-696 1 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

Office  of  the  Chief  of  Naval  Research 
800  North  Quincy  Street.  BCT  #  1 
Arlington,  VA  22217-5000 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

We  describe  analyses  and  numerical  computations  that  examine  the  validity  of  certain  common  simplifica¬ 
tions  in  the  chemical  and  thermophysical  model  on  the  calculation  of  detonation  properties.  In  particular,  we 
examine  the  effects  of  ( 1 )  approximations  of  the  temperature  dependence  of  the  specific  heat  C,  and  the  ratio  of 
specific  heats  7  =  C^/C,  in  a  detailed  numerical  simulation  of  detonation  in  a  mixture  of  hydrogen  and  oxygen 
mixture.  (2)  varying  the  heat  release,  and  (3)  varying  the  chemical  induction  time,  and  compare  all  of  those 
approximations  to  a  computation  that  uses  a  detailed  model  of  the  chemical  kinetics  and  correct  thermophysics. 
Finally,  we  describe  computations  of  detonations  in  mixtures  of  hydrogen  and  oxygen  diluted  with  argon  and 
carbon  tetrafluoride  and  show  how  the  numerical  model  correctly  predicts  the  trends  in  the  detonation  velocity 
as  the  percentage  of  carbon  tetrafluoride  is  increased. 


14.  SUBJECT  TERMS 

Detonation  structure 
Detonation  simulations 


Detonation  Modeling 


IS  NUMBER  OF  PAGES 
22 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  LIMITATION  OF  ABSTRACT 


NSN  7S40.01J80-5S00 


Slandatc)  form  298  (Rev  2-89) 

PtescriOM  bv  ANSI  Std  239  ’8 
298  102 


CONTENTS 


1.  Introduction  .  1 

2.  Reference  State  .  2 

3.  Calculation  of  Internal  Energy,  Pressure  and  Temperature  .  3 

3.1  Full  Integration  .  3 

3.2  Constant  Heal  Capacity  Approximation  .  3 

3.3  Comments  .  6 

4.  A  Detonation  in  a  Hi-Oi-Ar  Mixture  .  6 

4.1  Numerical  Model  .  6 

4.2  Cases  Studied  .  7 

4.3  Results  and  Comments  .  9 

5.  Detonations  in  Mixtures  Containing  Flurocarbons  .  10 

5. 1  Preliminary  Observations  .  10 

5.2  Time-Dependent  One-Dimensional  Calculations  .  11 

5.3  Comments  .  12 

6.  Conclusion  .  13 

Acknowledgments  .  13 

References  .  14 


iii 


Acossalon  Per 
MTIS 
DTir  T’.' 
Uun:'i70’  .1;  .  'id 
Ju.-tflc 


By...  _ _ _ 

[__Dl 

^  vail/ibl  1 1  i  y  .".ji'.os 
Ji\  V 1  i  tii.' I '  /  3 

JOlat,  i 


I 

H 


— I 


□□a 


Computations  of  Detonation  Structure: 
The  Influence  of  Model  Input  Parameters 


1.  Introduction 


Typical  reactive-flow  problems  solve  some  form  of  the  equations  of  conservation  of  mass,  momen¬ 
tum,  and  energy.  The  solution  of  these  equations  provides  the  knowledge  of  the  three  associated 
fundamental  quantities,  that  is,  density  p  (kg/m^),  momemtum  pv  (kg/m^s),  and  density  of  total 
energy  E  (J/m^).  The  total  energy  considered  is  usually  defined  as  the  sum  of  the  internal  energy 
and  the  kinetic  energy. 


(1) 


where  U  (J/kg)  is  the  internal  energy  of  the  system.  The  internal  energy  is  related  to  temperature 
and  pressure  through  the  caloric  equation  of  state. 


U  =  U(T,P),  (2) 

and  temperature  (K)  and  pressure  (Pa)  are  related  by  the  thermal  equation  of  state 

P  =  P{T,p).  (3) 

In  this  formulation,  variables  such  as  internal  energy,  temperature,  velocity,  and  pressure  are  derived 
quantities.  At  each  stage  of  the  calculation,  the  internal  energy  can  be  calculated  using  equation  (1), 
and  knowing  U  and  the  two  equations  of  state  enables  us  to  calculate  the  pressure  and  temperature. 
The  pressure  is  required  for  the  calculation  of  the  source  terms  in  the  conservation  equations.  For 
example,  the  temperature  is  required  when  chemical  and  thermal  processes  are  involved. 

The  calculation  of  the  energy  of  a  system  is  based  on  the  definition  of  a  ‘state’  at  which  the 
energy  is  assumed  to  have  a  given  value,  that  is,  a  ‘reference  state.’  Thus  we  first  define  this  concept 
as  accurately  as  possible  and  then  discuss  the  caloric  equation  of  state  in  its  simplest  form,  that  is, 

U  =  U{T)  =  Uf  +  I  C„[T)  dT  , 

where  Uj  \s  the  internal  energy  of  formation  and  the  integral  represents  the  sensible  internal  energy. 
This  form  leads  to  the  use  of  the  ratio  of  heat  capacities,  7  =  Cp/C„  and  automatically  limits  the 
use  of  the  material  in  this  report  to  mixtures  of  ideal  gases.  In  this  paper,  we  do  not  consider 
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other  more  complex  equatioDS,  such  as  U  =  U{p,T)  or  U{P,T).  Then,  the  relation  between  the 
temperature  and  the  pressure  is 

P(T,p)  =  pNRT, 

where  N  is  the  total  number  of  mol  (mol/kg). 

In  meiny  computations,  we  use  simplified  models  of  the  chemistry  for  which  we  do  not  track 
all  of  the  chemical  species  that  might  be  involved.  For  such  computations,  we  need  simplified 
models  of  Cv  The  question  is  then  how  to  get  the  best  model  or  one  that  is  consistent  with 
the  reduced  chemical  model.  The  purpose  of  this  report  is  to  describe  and  to  discuss  the  relations 
between  the  internal  energy  of  gaseous  mixtures  and  the  temperature  and  pressure,  that  is,  equation 
(2).  We  first  focus  on  theoretical  ways  to  compute  the  expression  J  C„{T)  dT  and  then  on  the 
temperature  dependences  of  the  heat  capacity,  C„,  and  the  ratio  of  heat  capacities,  7.  Then  the 
influence  of  the  approximation  to  Cv  on  the  derived  quantities  P  and  T  is  estimated.  In  the  second 
part  of  this  work,  the  influence  of  the  temperature  and  species  dependence  of  the  quantity  7  will 
be  studied  for  the  particular  case  of  a  one-dimensional  simulation  of  a  detonation  wave.  Finally, 
we  discuss  two  other  chemistry-dependent  parameters,  namely  the  heat-release  function  and  the 
chemical  induction  time. 

2.  Reference  State 

Energy  (enthalpy)  changes  are  relative  quantities,  thus  we  must  relate  each  compound  to  a  reference 
state  which  is  done  by  referring  each  energy  (enthalpy)  to  the  chemical  elements  in  their  standard 
reference  state.  By  definition,  we  assume  that  a  gaseous  mixture  is  in  the  reference  state  when; 

1.  All  constituents  are  in  the  standard  state,  that  is,  H2  for  hydrogen,  O2  for  oxygen,  N2  for 
nitrogen,  solid  carbon  graphite  for  carbon,  etc. 

2.  A  reference  temperature.  To,  is  chosen. 

3.  The  pressure  is  equal  to  1  bar  (100,000  Pa). 

In  this  reference  state,  the  internal  energy  (enthalpy)  of  the  mixture  is  assumed  to  be  equal  to 
zero.  We  may  thus  speak  about  the  ‘internal  energy  (enthalpy)  of  a  mixture’  instead  of  the  ‘change 
in  internal  energ>'  (enthalpy)  of  the  mixture.’ 
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3.  Calculation  of  the  Internal  Energy,  Temperature,  and  Pressure 
The  internal  energy  of  a  mixture  of  ideal  gases  is  given  by 

U{T)  =  +  r c.i{T)  dT  , 

i  i 

where  Cvi{T)  is  the  heat  capacity  of  the  species  i  at  constant  volume  (J/K  mol),  Ui  is  the  number 
of  moles  of  species  i  (mol /kg),  and  is  the  molar  internal  energy  of  formation  of  species  i  at 
temperature  To  (J/mol).  If  we  consider  the  energy  of  formation  of  the  whole  mixture,  C//,  and  its 
global  heat  capacity,  C„, 

f/;  =  and  C„(T)  =  ^ruCodT)  , 

1  t 

we  may  write 

U{T)  =  [/;+/  Co{T)dT  .  (4) 

JTo 

It  is  important  to  notice  that  Uj  is  the  energy  of  formation  of  the  mixture  at  the  chosen  reference 
temf>erature  Tg.  The  integral  in  equation  (4)  represents  the  sensible  internal  energy  of  the  system 
and  can  be  calculated  by  full  integration  or  by  assuming  a  temperature-independent  value  of  the 
heat  capacity  C„.  The  question  is  how  can  we  deduce  the  temperature  from  a  known  internal 
energy  U(T),  using  equation  (4). 

3.1  Full  Integration 

The  unknown  temperature  is  the  upper  boundary  of  the  integral  in  equation  (4).  Usually  C„(T)  is 
known  in  the  form  of  a  polynomial,  so  that  the  integral  can  be  evaluated  analytically  and  T  can  be 
calculated  iteratively.  In  this  case,  the  temperature  and  species  dependence  of  the  internal  energy 
is  fully  taken  into  account.  The  disadvantage  of  this  approach  is  the  time  it  may  take  to  iterate  to 
the  correct  solution  and  the  fact  that,  in  this  case,  the  densities  of  all  the  species  must  be  tracked. 
This  means  solving  many  complex  chemical  equations  and  increases  the  number  of  conservation 
equations. 

3.2  Constant  Heat  Capacity  Approximation 

In  order  to  save  computer  time,  it  is  often  convenient  to  use  simplified  chemical  models  and  a 
temperature-independent  heat  capacity  C„.  If  we  know  an  average  value  of  the  heat  capacity  of 
the  mixture  C*,  we  may  write 

U{T)  =  C:{T  -  To)  +  UI  ,  (5) 
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or  using  the  average  ratio  of  heat  capacities,  7*  =  C*ICl, 


RT, 


RT 


-  Wr^)  *  WF^)  ■ 


where  M  is  the  mean  molar  mass  of  the  gas  (g/mol).  The  total  energy  density  is  thus 

,,2 


(6) 


If  the  reference  temperature  is  equal  to  OK,  equation  (6)  can  be  simplified,  that  is, 


E 


=  pU°f  + 


7*  - 1 


(7) 


In  this  case  (To  =  0  K),  we  may  use  either  the  internal  energy  of  formation  or  the  enthalpy  of 
formation  Hj,  since  the  enthalpy  and  internal  energy  of  formation  are  equal  at  0  K  .  Equation  (7)  is 
particularly  simple  and  useful.  The  pressure  can  be  directly  calculated  from  equation  (7),  and  the 
temperature  from  the  thermal  equation  of  state.  The  theoretical  average  value  of  the  heat  capacity 
C*  and  its  derived  variable  7*  are  given  by 


c:  = 


C,{T)dT 


and 


7  = 


(8) 


If  this  value  of  7*  is  used,  the  pressure  computed  from  equation  (7)  must  be  done  using  the  full 
temperature  and  species  dependence  of  U{T).  As  U{T,ni)  is  usually  not  known  exactly  unless  the 
full  chemical  and  thermodynamic  computations  are  carried  out,  we  still  need  to  have  an  estimate 
of  C*  or  7*  for  simplified  chemical  models. 

We  denote  the  estimated  mean  values  of  C*  and  7*  by  and  y'*".  Then  the  value  of  7‘''  can 
be  estimated  theoretically 

a.  As  7+  =  1.66  for  a  monatomic  gas  and  'y'*'  =  1.40  for  a  diatomic  gas, 

b.  By  calculating  'y'^  at  room  temperature  (298  K)  and  assuming  that  7+  is  constant,  so  that 


_4-  E.«iC„i.298  + 

7  =  - ^=1 -  =  7298 

2^i  n«Cri.298 


(9) 


c.  By  assuming 


_  C„T  +  .  +  _  C+  -t-  R/M 


1 


ilT  -  1)  ’  +  (7rM8  “  I)  ’ 


(10) 


(11) 


(7+  -  1)  2 

if  a  reasonable  value,  Cvt  or  7t,  of  Cv  or  7  at  the  unknown  temperature  T  can  be  estimated. 
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For  the  species  O2  and  H2O,  Figures  1  and  2  show  four  quantities  computed  using  the  different 
approximations  of 

A.  the  value  of  at  298  K, 

B.  C„  calculated  at  the  current  temperature  T,  that  is,  CvT, 

C.  the  rigorous  value  C*  calculated  with  equation  (8),  and 

D.  the  mean  value  of  calculated  with  equation  (10). 

The  first  graphs,  Figure  la  or  2a,  show  the  heat  capacity  calculated  at  constant  volume  Cy  The 
second  graphs,  Figure  la  or  2b,  show  the  ratio  of  heat  capacities  computed  using  the  same  approx¬ 
imations.  All  calculations  assume  To  —  OK. 

In  addition  to  the  standard  quantities  and  7,  Figures  1  and  2  also  show  two  other  quantities, 
errorl  and  error2.  From  equation  (7),  we  calculate  the  error  in  the  pressure  when  the  value  7'*'  is 
used  instead  of  the  exact  7* .  This  relative  error  is  given  by  the  ratio 

errorl  ==  — - ^  .  (12) 

7*  —  1 

The  quantity  errorl  also  represents  the  relative  error  in  the  calculated  temperature  because  we  are 
assuming  ideal  gases  and  the  ideal  gas  law.  The  curves  shown  in  Figures  Ic  and  2c  are  calculated 
with 


A.  7''‘  =  7298,  as  in  equation  (9), 

B.  7"^  =  77>,  that  is,  7(T)  at  the  current  temperature,  and 

C.  7''‘  from  equation  (11). 

Equation  (11)  might  appear  to  be  a  good  estimation  of  7*,  but  to  compute  7"^  using  this  equa¬ 
tion  requires  knowledge  of  the  temnerature  in  the  term  7r.  In  order  to  determine  how  accurately 
7x  must  be  known  in  equation  (11),  we  define  the  quantity  error!  in  the  following  way.  At  the 
current  temperature  Tc,  we  calculated  errorl  from  equation  (12)  by  using  7+  from  equation  (11), 
but  now  for  three  different  temperatures,  T,  and  T  ±  200.  We  plot  the  maximum  of  these  three 
values 


error2{Tc,  ')Tc  ) 


max 


Vi+m-n 

/7+(r)-n 

h^{T)-\\ 

\  7*  ~  1  /  T=Tc-200 

\  7*  ■“  1  /  T-Tc 

V  7*  ~  I  /  T=Tc+200 

Figures  Id  and  2d  thus  show  the  largest  error  in  pressure  calculated  with  equation  (7),  when  the 
temperature  required  to  calculate  the  quantity  77-  in  equation  (11)  is  known  w'ithin  a  range  of 
±  200  K.  The  curves  labelled  ‘±  500’  in  the  error!  graphs  have  the  same  definition  for  range  of 
temperatures  ±  500. 
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3.3  Comments 


The  error  in  the  pressure  or  temperature  calculated  with  equation  (7)  and  using  a  constant  value  of 
7  can  be  large.  When  calculations  are  carried  out  with  7298,  the  error  goes  up  to  25%  for  diatomic 
species  and  up  to  50%  for  H2O.  The  inaccuracy  becomes  especially  significant  for  temperatures 
higher  than  1000  K.  The  same  error  using  the  mean  value  7'*'  from  equation  (11)  is  less  than  5% 
in  the  whole  range  of  temperature  above  1000  K. 

The  use  of  an  estimated  temperature-dependent  7"^  instead  of  a  constant  7293  is  more  important 
when  tri-atomic  or  multi-atomic  molecules  are  involved  (see  H2O)  because  when  the  molecule  has 
more  degrees  of  freedom,  the  heat  capacity  has  a  stronger  temperature  dependence.  For  nonreactive 
flow,  the  assumption  of  constant  7298  is  often  satisfactory.  Usually  this  kind  of  calculation  involves 
air,  that  is  a  mixture  of  diatomic  species,  and  the  temperature  does  not  increase  as  much  as  in  a 
reactive  flow. 

When  the  phenomena  studied  are  highly  temperature  dependent  (for  example,  exothermic 
chemically  reacting  flow),  an  estimation  of  7*  is  necessary.  Depending  on  the  case  studied,  a  rough 
estimation  of  Jt  and  the  use  of  equation  (11)  produces  satisfactory  estimates. 

The  evaluation  of  the  quantity  errorl  shows  that  the  calculation  of  pressure  using  equation  (7) 
is  very  sensitive  to  the  temperature  dependence  of  7.  Not  including  the  temp)erature  dependence 
can  lead  to  an  extremely  inaccurate  estimate  of  7.  However,  for  practical  applications  in  the 
evaluation  of  pressure,  the  temperature  dependence  of  7r  does  not  have  to  be  very  accurate;  a 
rough  estimate  is  quite  satisfactory,  as  shown  from  the  values  of  error2. 

4.  A  Detonation  in  a  H2-02-Ar  Mixture 

In  this  section  we  will  study  the  effect  of  including  the  temperature  dependence  of  7  and  the 
influence  of  the  heat-release  function  on  a  computed  detonation  wave. 

4.1  Numerical  Model 

The  model  solves  the  three  basic  conservation  equations  (mass,  momentum,  and  energy)  with  a 
model  for  chemical  reactions  and  heat  release.  The  chemical  model  is  a  two-step  model  where  the 
steps  correspond  to  the  chemical  induction  period  and  the  heat-release  period. 

Chemical  Induction  Period 

We  keep  track  of  the  evolution  of  the  chemical  induction  time  by  solving  an  equation  of  the 

form 

dTind  _  1 

“IT  "  uUT{t),P{t))  ' 
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where  Tind  represents  the  fraction  of  elapsed  induction  time  and  0<  Ti„d  <  1.  The  induction  time 
tind  is  calculated  for  the  local  conditions  of  temperature  and  pressure  and  may  be  obtained  by 
integrating  a  detailed  set  of  reaction  rates  or  from  experimental  data.  In  our  case,  we  have  used  a 
set  of  data  provided  by  Oran  et  al.^  During  the  induction  period,  no  energy  is  released. 

Heat-Release  Period 

The  exothermic  chemical  reaction  is  tracked  by 

d^hrel  1 

dt  threl  ’ 

where  Threi  represents  the  fraction  of  released  energy  and  0<  Threi  <  1-  This  progress  parameter 
Threl  is  thus  equal  to  AU^{t)/AU%  The  ‘heat-release  time’  (threi)  is  either  kept  constant  over 
the  entire  heat-release  period  or  may  evolve  as  a  function  of  Threi-  The  function  threii'^hrei)  may 
be  obtained  from  kinetic  calculations  that  integrate  a  set  of  elementary  chemical  reactions  subject 
to  the  constraints  of  steady-state  shock  hydrodynamics. 

4.2  Cases  Studied 

The  specific  mixture  studied  is  a  H2:02:Ar  /  2:1:7,  at  an  initial  pressure  and  temperature  of  50  Torr 
and  298  K,  respectively.  We  have  performed  calculations  for  the  eight  different  cases  summarized 
in  Table  1.  Note  that: 

1.  The  heat-release  time  threi  is  implemented  in  three  different  ways,  that  is, 

a.  fast  reaction  with  constant  rate:  threi  =  5  ns, 

b.  slow  reaction  with  constant  rate:  threi  =  20  ns, 

c.  variable  rate  threi{Threi)  provided  by  a  chemical  kinetic  computation  in  which  the  total 
heat-release  time  is  15-17  /xs. 

2.  The  maximum  amount  of  energy  released  is  provided  by  the  chemical  kinetic  computation, 
except  in  case  6.  This  case  is  discussed  later. 

3.  The  heat  capacity  ratio  7  is  either  constant,  in  which  case  7  is  equal  to  its  initial  value  of  1.556, 
or  7  evolves  as  a  function  of  the  progress  parameter  t,  taking  into  account  a  temperature  and 
species  dependence  of  the  heat  capacity.  This  dependence  is  estimated  by  a  preliminary  kinetic 
computation. 

4.  The  mean  molecular  mass,  M,  is  either  constant  or  a  function  of  r,  as  is  7. 

5.  The  induction  time  is  tabulated  as  a  function  of  temperature  and  pressure  by  integrating  a  full 
set  of  reaction  rates.  The  purpose  of  cases  7  and  8  is  to  show  the  influence  of  the  induction 
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Table  1.  Physico-Chemical  Parameters  Used  in  Simulations 


Case 

^krel 

7 

M 

^xnd 

[ms] 

IkJ/kg] 

(1) 

(2) 

(3) 

(4) 

(5) 

1 

const  (5  /is) 

1010 

const 

const 

- 

2 

const  (20  /is) 

1010 

const 

const 

- 

3 

K'^hrel) 

1010 

const 

const 

- 

4 

const  (20  fis) 

1010 

f('^/ire/) 

f(^/irez) 

- 

5 

^{Threl') 

1010 

f(TZirc/) 

- 

6 

const  (10  fis) 

820 

const 

const 

- 

7 

fC'Oirc/) 

1010 

f(Thre/) 

f(Threz) 

X  1.5 

8 

f(T?ircz) 

1010 

f(^/irez) 

X  2.0 

9 

the  fluid  dynamics  equations  are  solved  with 

a  full  set  of  reaction  rates 

and  self-consistent  heat  capacities. 

time  on  the  detonation  propagation,  so  we  have  used  the  same  conditions  as  in  case  5  but 
multiplied  the  tabulated  induction  time  by  factors  of  1.5  and  2.0,  respectively. 

6.  Case  6  describes  a  calculation  where  all  parameters  {threh  7,  and  M)  are  kept  constant  and 
the  total  amount  of  energy  released  is  artifically  reduced. 

For  each  case,  we  focus  on  the  detonation  velocity,  the  thickness  of  the  induction  zone  and  the 
thickness  of  the  total  reaction  zone.  The  detonation  velocity  may  be  compared  to  the  Chapman- 
Jouguet  (CJ)  value,  that  is,  1618  m/s. 

In  summary, 

-  Cases  1,  2,  and  3  keep  7  and  M  constant,  and  show  the  influence  of  heat  release. 

-  Cases  4  and  5  show  the  influence  of  the  temperature  and  species  dependence  of  7  and  A/. 

-  Cases  7  and  8  show  the  influence  of  the  induction  time. 

-  Case  9  is  a  complete  reference  calculation  with  a  full  solution  of  the  coupled  fluid-chemical 
problem  and  a  detailed  set  of  elementary  reaction  rates. 
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4.3  Results  and  Comments 


The  initial  conditions  of  the  problem  are  a  typical  shock-diaphragm  problem  where  the  low  pressure 
is  50  Torr  (6666  Pa)  and  the  driving  gas  is  at  an  initial  pressure  and  temperature  of  2280  Torr 
(304,000  Pa)  and  2300  K,  respectively.  After  a  few  time  steps,  the  mixture  ignites  and  the  shock 
velocity  increases  suddenly  and  strongly.  After  an  initial  relaxation,  the  velocity  stabilizes  and  the 
result  is  a  coupled  shock  front  and  reaction  zone.  Table  2  gives  the  values  of  the  steady  detonation 
velocity,  D,  the  thickness  of  the  induction  zone,  Xi„a,  and  the  reaction  zone  Xreact  for  each  case 
described  above.  Figures  3  shows  results  as  a  function  of  time  for  case  5.  The  maximum  energy 
released  in  the  calculation  carried  out  with  the  full  chemistry  (case  9)  is  about  1080  kJ/kg,  that  is, 
a  little  more  than  the  value  used  in  the  other  calculations. 


Table  2.  Detonation  Velocity  and  Structure  of  Simulations 


rase 

1 

2 

3 

4 

5 

6 

7 

8 

9 

D  (m/s) 

1802 

1787 

1798 

1649 

1650 

1626 

1652 

1655 

1G70 

Xind  (mm) 

0.71 

0.77 

0.75 

1.2 

1.2 

1.2 

1.8 

2.2 

0.9/1.5 

Xreact  (mm) 

4.1 

13.3 

12.2 

12.9 

11.8 

7.4 

11.8 

12.5 

- 

Despite  a  significant  difference  in  the  heat-release  time,  the  computations  do  not  show  any 
appreciable  differences  on  the  detonation  velocity  D  among  cases  1,2,  and  3  and  between  cases  4 
and  5.  However,  D  is  quite  high  when  the  heat  capacity  ratio  is  kept  constant.  The  induction  zone 
is  not  affected  by  the  heat  release  function  and  the  total  reaction-zone  thickness  evolves  consistently 

with  threl- 

The  temperature  and  species  dependence  of  7  and  M  modifies  D  dramatically  (cases  4  and 
5  compared  to  1,  2,  and  3).  It  is  important  to  notice  that  allowing  7  and  M  to  be  variable 
produces  velocities  close  to  the  Capman-Jouguet  detonation  velocity  without  any  modification  of 
the  estimated  amount  of  released  energy.  The  results  obtained  with  7  =  /(r,  5  )  agree  well  with 
the  reference  calculation,  case  9. 

Case  6  shows  that  it  is  possible  to  obtain  the  same  qualitative  result  by  keeping  7  constant 
and  decreasing  the  amount  of  released  energy  (in  fact,  by  decreasing  the  fraction  of  reacted  fuel). 
A  reduction  of  the  total  heat  released  compensates  for  the  thermal  ideality  of  the  mixture. 

A  relatively  large  variation  of  the  induction  time  does  not  affect  the  steady  detonation  velocity. 
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However,  the  thickness  of  the  induction  zone  is  modified  consistently  by  the  change  in  induction 
time  (cases  7  and  8). 

The  computations  presented  above  and  the  results  obtained  are  valid  for  this  particular  hydro¬ 
gen-oxygen-argon  mixture,  and  the  observations  and  conclusions  could  be  reasonably  extended  to 
similar  mixtures.  We  are  not  sure,  however,  how  the  results  generalize  to  other  reactive  systems. 
To  help  answer  to  this  question,  we  now  describe  calculations  for  an  “exotic”  mixture  for  which 
experimental  data  are  available. 

5.  Detonations  in  Mixtures  Containing  Flurocarbons 

The  mixture  considered  is  stoichiometric  H2-O2,  diluted  50%  with  Ar  and  CF4,  at  an  initial  pressure 
and  temperature  of  200  Torr  (26400  Pal  and  298  K,  respectively.  The  computations  have  been 
carried  out  with  0,  5,  10,  15,  and  20%  of  CF4.  Here  there  are  many  more  kinds  of  chemical 
reactants  and  products  than  in  the  simple  mixture  of  H2  and  O2  discussed  above.  A  possible 
kinetic  mechanism  and  reaction  rate  constants  are  described  by  Nzeyimana  and  Van  Tiggelen^.  The 
advantage  of  testing  our  detonation  computation  using  this  mixture  is  that  we  have  an  extended  set 
of  experimental  data  provided  by  the  Laboratoire  de  Physico-Chimie  de  la  Compustion,  Louvain- 
La-Neuve,  Belgium.  Moreover,  an  interesting  feature  of  fluorocarbon  compounds  is  that  their  heat 
capacity  is  much  more  temperature  dependent  than  all  of  the  species  involved  in  the  pure  H2-O2 
mechanism. 


Table  3.  Experimental  and  CJ  Detonation  Velocities 

Composition  of  Mixture:  33.33%  H2,  16.66%  O2,  (50.00  -  A')%  Ar,  A''%  CF4 


X  (%  CF4) 

0.0 

5.0 

10.0 

15.0 

20.0 

Dcj  (m/s) 

1874 

1875 

1874 

1872 

1693 

Dcxp  (m/s) 

1830 

1859 

1854 

1829 

1760 

5.1  Preliminary  Observations 

The  behavior  of  explosive  mixtures  containing  additional  fluorocarbons  is  of  interest  because  these 
compounds  are  thought  to  inhibit  exothermicity  in  chemical  reactions.  However,  there  are  some 
that  behave  more  like  promotors  of  detonation,  such  as  CF4,  when  less  than  15%  is  added  to  the 
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reference  mixture  H2:02:Ar  /  2:1:3.  Moreover,  there  is  a  notable  difference  betwo^n  the  experimen¬ 
tal  detonation  velocity  and  computed  CJ  detonation  velocity.  Table  3  summarizes  some  of  these 
results. 

The  major  comment  about  the  experimental  results  shown  in  Table  3  is  that  additional  CF4 
slightly  increases  the  detonation  velocity,  an  effect  that  is  not  predicted  by  CJ  calculations.  More¬ 
over,  when  a  large  amount  (20%)  of  CF4  is  added,  the  experimental  detonation  velocity  becomes 
surprisingly  larger  than  the  expected  CJ  value.  To  better  understand  the  CF4  behavior  in  detona¬ 
tion  waves,  calculations  were  done  to  test  the  chemical  mechanism.^  It  is  important  to  note  that  the 
kinetic  mechanism  of  a  reactive  system  involving  CF4  is  not  well  known  and  most  of  the  reaction 
rates  involving  CF4  or  one  of  its  decomposition  products  are  estimated.  Nevertheless,  these  calcu¬ 
lations  showed  that  the  products  of  the  chemical  reaction  tracked  by  a  kinetic  mechanism  are  quite 
different  than  those  predicted  by  a'  y  equilibrium  calculations.  One  of  the  main  assumptions  of  the 
CJ  calculation  is  that  the  chemicf  reactions  are  in  an  equilibrium  state,  and  it  is  this  assumption 
that  could  explain  why  the  CJ  theory  fails  in  this  particular  case.  Taking  this  shift  in  composition 
into  account,  a  larger  amount  of  energy  is  released  in  the  nonequilibrium  kinetic  calculation  than 
in  a  calculation  based  on  an  equilibrium  compKisition. 


Table  4.  Computed  Detonation  Velocity  and  Structure  of  Reaction  Zone 


Con.position  of  Mixture: 

33.33%  H2, : 

16.66%  O2,  (50.00 

-  X)%  At,  X%  CF4 

X  (%  CF4) 

0.0 

5.0 

10.0 

15.0 

20.0 

Dcalc  (m/s) 

1875.1 

1923.2 

1918.8 

1871.5 

1717.9 

^ind  (rnm) 

0.53 

0.48 

0.45 

0.48 

0.52 

^krei  (nirn) 

1.09 

1.08 

1.06 

1.19 

1.13 

5.2  Time-Dependent  One-Dimensional  Calculations 

The  detonation  calculations  are  performed  with  the  numerical  model  described  above  that  takes 
account  of  the  evolution  of  the  ratio  of  the  heat  capacities,  7,  and  the  mean  molar  mass,  M  as  a 
function  of  the  reaction  progress  parameter  Threi-  The  he-it  release,  7,  and  M  are  related  to  Threi 
by  a  set  of  tabulated  data  provided  by  kinetic  computation.  In  the  mixture  containing  up  to  15% 
CF4,  the  initial  pressure  and  temperature  in  the  driving  gas  section  are  904,000  Pa  and  3000  K, 
respectively.  In  the  mixture  containing  20%  CF4,  the  calculation  is  performed  with  pressure  and 
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temp)erature  equal  to  1,204,000  Pa  and  3300  K,  respectively.  Table  4  summarizes  the  results.  Figure 
4  shows  results  of  selected  computations  in  more  detail,  and  Figure  5  summarizes  the  detonation 
velocities  predicted  for  different  initial  concentrations  of  CF4. 

5.3  Comments 

The  absolute  differences  between  calculated  detonation  velocities  and  the  experimental  values  may 
be  due  to  a  number  of  factors,  such  as  the  geometry  of  the  shock  tube  or  multidimensional  effects. 
However,  more  significant  is  the  way  both  the  computed  and  experimental  detonation  velocities 
evolve  with  the  addition  of  CF4.  With  up  to  15%  CF4,  there  is  good  agreement  between  the 
experimental  data  and  the  calculation.  Moreover,  the  calculations  help  us  to  better  understand  the 
effect  of  adding  CF4. 

The  increase  of  the  detonation  velocity  is  the  consequence  of  two  competing  effects: 

1.  Addition  of  CF4  causes  an  increase  of  the  specific  energy  released  compared  with  the  reference 
mixture  H2:02:Ar. 

2.  The  temperature  behind  the  leading  shock  (von  Neumann  temperature)  decreases  when  CF4 
is  added  because  the  heat  capacity  of  mixtures  containing  CF4  and  its  decomposition  products 
(CF3,  CF2O,  ...)  is  large. 

The  chemical  model  enables  us  to  describe  the  interplay  of  these  two  effects  as  the  amount  of  CF4 
varies.  This  is  particularly  interesting  because  the  model  does  not  require  the  evolution  of  all  par¬ 
ticipating  chemical  compounds  to  be  followed.  A  elementary  kinetic  mechanism  of  reactive  systems 
like  H2-O2-CF4  involves  at  least  14  species  and  the  integration  of  such  a  chemical  mechanism  is 
demanding  in  computer  time.  The  trend  observed  when  the  mixture  contains  20%  is  less  satisfac¬ 
tory.  One  reason  for  this  may  be  our  poor  knowledge  of  the  chemistry  of  mixtures  containing  CF4 
and  the  fact  that  the  simulation  is  very  sensitive  to  the  total  heat  release  and  to  the  variation  of 
7,  as  shown  in  the  first  section  of  this  work. 

Mixtures  with  large  amounts  of  CF4  are  also  more  difficult  to  ignite,  this  is  one  of  the  effects 
of  the  large  heat  capacity  of  CF4.  The  structure  in  the  reaction  zone  fluctuates  noticeably  when 
there  is  over  10%  of  CF4.  Whether  these  oscillations  are  from  a  numerical  origin  or  are  a  physical 
property  of  the  system  has  not  yet  been  investigated.  However,  they  may  be  related  to  the  increasing 
irregularity  observed  in  multidimensional  cellular  structure  of  mixtures  containing  over  10%  CF4. 
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6.  Conclusion 


This  work  discusses  the  accuracy  of  selected  input  data  used  in  computations  of  detonations. 
In  particular  we  have  considered  the  effect  of  including  the  temperature  dependence  of  the  heat 
capacity  in  the  numerical  model.  The  magnitude  of  the  error  made  with  constant  7  is  particularly 
large  when  the  temperature  is  higher  than  1000  K.  This  error  increases  substantially  when  chemical 
reactions  modify  the  composition  of  the  mixtures  studied  and  produce  multi-atomic  compounds. 
A  modified  chemical  model  has  been  used  in  order  to  include  a  temperature-dependent  7  in  the 
computation.  The  model  consists  of  a  two-step  chemical  process  but  takes  into  account  a  relatively 
detailed  description  of  the  elementary  chemical  reaction  rates.  The  model  has  been  tested  for 
the  mixture  H2:02:Ar  and  compared  successfully  with  other  chemical  models.  Calculations  of 
detonations  in  mixtures  containing  more  complex  species  and  involving  a  more  complex  chemical 
mechanism  have  been  tested  and  show  a  good  qualitative  agreement  with  available  exp)erimental 
data.  We  therefore  believe  that  this  modified  chemical  model  can  be  used  for  other  reactive  systems 
whose  chemical  kinetics  is  known  to  within  a  reasonable  accuracy. 

Some  points  remain  unanswered.  We  do  not  know  to  what  extent  the  kinetic  calculations 
used  to  fit  the  set  of  data  dr/dt,  7,  and  M  may  influence  the  fluid  dynamic  model.  We  have 
used  a  kinetic  mechanism  related  to  the  Raleigh  line  of  the  steady  shock  theory,  but  it  would  be 
interesting  to  investigate  the  effect  of  isothermal  kinetics  or  kinetics  at  constant  density.  Two- 
dimensional  calculations  should  be  a  good  test  of  how  accuratly  this  model  describes  the  chemistry 
occurring  behind  the  shock  wave. 
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Figure  1.  a)  Heat  capacity  at  constant  volume  (J/mol  K),  b)  ratio  of  specific  heats,  7.  c)  errorl 
(defined  in  text)  and  d)  error2  (defined  in  text)  as  a  function  of  temperature  for  Oa- 
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Temperacure  (K)  Temperature  (K) 


Figure  2.  a)  Heat  capacity  at  constant  volume  (J/mol  K),  b)  ratio  of  specific  heats,  7,  c)  error  1 
and  d)  error'!  as  a  function  of  temperature  for  H2O. 


Figure  3.  a)  Detonation  velocity  and  b)  reaction  zone  as  a  function  of  time  derived  from  a  one¬ 
dimensional  simulation  of  a  detonation  propagating  in  a  mixture  of  H2:02:Ar/2:l:7  at  50  Torr  and 
298  K. 
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Figure  4.  a)  Detonation  velocity  and  b)  reaction  zone  as  a  function  of  time  derived  from  a  one¬ 
dimensional  simulation  of  a  detonation  propagating  in  a  H2-O2  mixture  containing  10%  CF4  at 
200  Torr  and  298  K. 


18 


Figure  5.  Numerical,  experimental,  and  CJ  detonation  velocities  in  a  H2-O2  mixture  as  a  function 
of  the  amount  of  added  fluorocarbon. 
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